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Abstract: Retrievals of ice cloud properties using infrared measurements at 3.7, 6.7, 7.3, 
8.5, 10.8, and 12.0 pm can provide consistent results regardless of solar illumination, but 
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are limited to cloud optical thicknesses x < ~6. This paper investigates the variations in 
radiances at these wavelengths over a deep convective cloud system for their potential to 
extend retrievals of x and ice particle size D e to optically thick clouds. Measurements from 
an imager, an interferometer, the Cloud Physics Lidar (CPL), and the Cloud Radar System 
(CRS) aboard the NASA ER-2 aircraft during the NASA TC 1 * * 4 (Tropical Composition, 
Cloud and Climate Coupling) experiment flight during 5 August 2007, are used to examine 
the retrieval potential of infrared radiances over optically thick ice clouds. Simulations 
based on coincident in situ measurements and combined cloud x from CRS and CPL 
measurements are comparable to the observations. They reveal that brightness 
temperatures at these bands and their differences (BTD) are sensitive to x up to ~20 and 
that for ice clouds having x > 20, the 3.7-10.8 pm and 3. 7-6. 7 pm BTDs are the most 
sensitive to D e . Satellite imagery appears to be consistent with these results suggesting that 
x and D e could be retrieved for greater optical thicknesses than previously assumed. But, 
because of sensitivity of the BTDs to uncertainties in the atmospheric profiles of 
temperature, humidity, and ice water content, and sensor noise, exploiting the small BTD 
signals in retrieval algorithms will be very challenging. 

Keywords: clouds; optical depth; particle size; satellite; TC 4 ; multispectral thermal infrared 


1. Introduction 

Because of their influences on the Earth’s radiation budget and hydrological cycle, the need to 

monitor cloud properties from satellites has long been recognized in the study of climate (e.g., [1-5]). 
With the widespread availability of well-calibrated polar and geostationary orbiting satellite imagers, 
more sophisticated cloud parameterizations, maturing remote sensing algorithms, and faster computers, 
cloud properties sensed remotely from satellite imager data are also increasingly relevant for improving 
numerical weather forecasts (e.g., [6,7]). Because of their favorable spatial and temporal coverage, 
passive satellite measurements at solar and infrared (IR) measurements have been extensively used to 
retrieve optical and micro/macro physical properties of clouds. A wide variety of techniques have been 
developed to retrieve cloud properties using various combinations of those radiances measured by both 
operational and research satellite imagers such as the Moderate Resolution Imaging Spectroradiometer 
(MODIS). Most methods have focused on retrieving cloud properties, such as optical thickness x and 
particle effective diameter D e , during the daytime because reflected solar radiation can be used to 
retrieve the properties over nearly the full range of expected values (e.g., [4,8-13]). During both day 
and night, it is possible to estimate cloud top heights for the full range of cloud optical thicknesses 

(e.g., [5,14,15]), but retrievals of x itself and D e have been limited to optically thin (x < ~5) clouds 
because of the constraints of the blackbody limit (e.g., [10,16-24]). Thus, monitoring of cloud x and 
D e , parameters that determine the cloud water path CWP, using satellite solar and IR imagers is 
confined to daytime for the clouds that contribute the bulk of liquid and ice water to the atmospheric 
water budget. 
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The diurnal variations of both liquid and ice cloud properties are also important for understanding 
the Earth radiation and latent heat budgets. Accurately quantifying the diurnal variations of cloud 
properties then is crucial for improving climate models (e.g., [25-29]). Thus, retrieval of a full range 
of cloud properties during nighttime is essential for accurately defining their role in climate. Moreover, 
without that full range, assimilations of x or CWP will produce diminished benefit to numerical 
weather predictions at night compared to those during the day. Mitigation of the blackbody limitation 
in cloud remote sensing, therefore, is an endeavor important for further validation and improvement of 
both weather and climate models. 

To begin the process of addressing this challenge, this study investigates the IR signals from 
optically thick ice clouds using campaign measurements and rigorous radiative transfer modeling to 
determine the potential for estimating optically thick ice cloud x and D e using selected IR bands at 3.7, 
6.7, 7.3, 8.5, 10.8, and 12.0 pm. These bands are among the channels on many geostationary and 
polar-orbiting meteorological satellite imagers. Section 2 describes the detailed measurements from the 
NASA ER-2 and DC-8 aircraft and radiosonde ascents for a case study of a tropical deep convective 
cloud system observed during the Tropical Composition, Cloud and Climate Coupling Experiment 
(TC 4 ) [30], The measurements include cloud vertical profiles made by a cloud radar and lidar, ice 
cloud particle sizes and habits from in situ probes, and radiances made by an airborne scanning 
spectrometer covering visible to IR bands and an airborne scanning high-resolution IR interferometer. 
Using observations and simulations, section 3 examines the effects of x and D e on upwelling IR 
radiances over the tropical deep convective cloud system. The measurements from the airborne scanning 
high-resolution IR interferometer are used to construct IR brightness temperatures at the MODIS bands 
of 3.7, 6.7, 7.3, 8.5, 10.8, and 12.0 pm. Simulations are based on a rigorous radiative transfer model with 
inputs of in situ atmospheric profiles and cloud microphysical and optical properties derived from 
airborne measurements. Gamma particle size distributions are also used in simulations for the aim of 
investigating the influence of D e on IR radiances. Section 4 summarizes this study. 

2. Combined Measurements over a Deep Convective Cloud System 

This study focuses on a particular convective cloud system observed during the 5 August 2007 TC 4 
mission, when the ER-2 and DC-8 executed more than 4 hours of flights at their cruising altitudes of 
about 20 and 10 km, respectively. The coordinated flight segment between 1350 and 1438 UTC is the 
focus of this study. 

Because of its high altitude, the ER-2 effectively serves as a satellite having an array of remote 
sensing instruments. The MODIS/Advanced Spacebome Thermal Emission and Reflection Radiometer, 
ASTER Airborne Simulator (MASTER) on the ER-2 measured radiances at most of the satellite 
imager wavelengths of interest. The MASTER is a cross-track scanning spectrometer that measures 
reflected solar and emitted thermal radiation in 50 narrowband channels in the range 0.4-13 pm. Its 
2.5-mrad instantaneous field of view (FOV) results in a spatial resolution at nadir of 50 m at the 
nominal ER-2 altitude of 20 km. At that altitude, the total FOV of the MASTER is 85.92° (±42.96° 
from nadir), which yields a ground swath of about 37.25 km. The MASTER is described in detail 
by [31] and [32], The IR radiances observed by MASTER bands 30 (3.7 pm), 43 (8.6 pm), 48 
(1 1.3 pm), and 49 (12.2 pm) are used in this study. 
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During TC 4 the MASTER lacked the water vapor bands at 6.7 and 7.3 pm that are currently carried 
on the MODIS. In this study, the measurements from the Scanning High-resolution Interferometer 
Sounder (S-HIS) on the ER-2 are used to construct radiances in the 6.7 and 7.3 pm bands that are 
available on the MODIS. The S-HIS is a scanning interferometer which measures emitted thermal 
radiation at high spectral resolution between 3.3 and 18 pm [33]. Its instrument FOV is 100 mrad and 
its total view angle is ±30°, so from the ER-2 it produces sounding data with 2 km resolution at the 
nadir and covers a 40-km ground swath. 

The ER-2 carried the CRS and CPL, which are well suited for ice cloud studies because of their 
high sensitivities and spatial resolutions. The CRS [34] is a 94-GHz (W-band, 3.2-mm wavelength) 
pulsed polarimetric Doppler radar. It was designed to provide nadir cloud profiling from the ER-2. The 
ER-2 CRS measurements have a vertical resolution of 37.5 m and temporal resolution of 0.5 s. The 
minimum detectable reflectivity is -28 dBZ for the CRS at a range of 15 km from the ER-2. Thus, the 
CRS measurements could miss the tops of ice clouds with small ice particles. Lidar measurements are 
used to observe the clouds missed by the CRS since lidar is more sensitive to small ice particles and 
optically thin cirrus than radar. McGill et al. [35] found that the x missed by the CRS, but detected by 
the CPL, is in the range of 0.15-0.45. The CPL is a backscatter lidar with three wavelengths of 0.355, 
0.532, and 1.064 pm [36]. It scans downward from the ER-2 with a vertical resolution of 30 m and a 
1-s temporal resolution, corresponding to -200 m at an average ER-2 ground speed of 200 ms -1 . Only 
the CPL measurements at 0.532 pm over the deep convective cloud system are used in this study. 

The DC-8 carried several cloud particle size probes to take in situ measurements and dropsondes to 
obtain temperature and humidity profiles below the aircraft. Ice cloud bulk scattering properties are 
mainly dependent on ice particle habits and ice cloud particle size distributions. Combined 
measurements from the 2-Dimensional Stereo (2D-S) probe, 2D Cloud Imaging Probe (CIP), and 
Precipitation Imaging Probe (PIP) on the DC-8 are used to construct ice cloud particle size 
distributions covering size bins from 5 pm to 6 mm in the present study. 

Numerous small ice particles (maximum dimension D less than 50 pm) have been observed in ice 
clouds by in situ measurements. However, large portions of these small ice particles could be 
generated by large crystals shattering on the inlets of instruments such as the Forward Scattering 
Spectrometer Probe (FSSP) and others which ingest ice particles into sample volumes through inlets 
([37] and references therein). The 2D-S uses detector arms instead of an inlet tube to minimize 
shattering fragments from entering the sample volume, however, collisions of the ice crystals and the 
forward elements of the probe still introduce some significant uncertainties into the derived size 
distributions [38,39], It detects particles as small as 10 pm (the minimum size bin is 5-15 pm) [38], 
For the TC 4 dataset, the maximum particle size bin that produced adequate sampling statistics was 
about 3000 pm. 

The CIP measured particle sizes ranging from about 50-100 pm to more than 1 mm, and the PIP 
measured particles sizes from about 100 pm to 6 mm [40]. Although both instruments nominally 
provide a minimum detectable size of 50 pm, the lowest detection threshold from the CIP is 100 pm 
during TC 4 because of uncertainties in the probe’s sample volume [40-42], Jensen et al. found that the 
size distributions derived from the CIP agree well with those from the 2D-S imaging instruments for 
the size bins above 100-200 pm despite the uncertainties in both instruments’ sample volume 
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sizes [37]. For detailed information on the CIP and PIP, see Tian et al. who investigated cirrus ice 
particle size distributions using TC 4 CIP and PIP observations [42], 

The Cloud Particle Imager (CPI) can image particles in the size range of 15-2,500 pm with a nominal 
2.3-pm resolution [38], During TC 4 , the CPI aboard the DC-8 had 33 size bins between 0 and 1,000 pm. 
The particle densities of seven habit categories including sphere, column, plate, rosette, budding 
rosette, small and large irregular particles were measured separately for ice particles with D < 50 pm 
and D > 50 pm. The habit mixtures derived with an automated technique from the measured number 
densities are used for computing ice cloud bulk scattering properties. Although all of the probe 
measurements have some significant uncertainties, they still provide the most realistic estimate of the 
microphysical properties of the clouds in this system and should yield reasonable radiances. 

Figure 1 . Infrared (10.8 pm) brightness temperatures from combined GOES-10 and 
GOES-12 data over the TC 4 domain, 5 August 2007. (a) 1558 UTC, with flight tracks for the 
day’s entire mission, ER-2 (blue) and DC-8 (red). Flight altitudes are shown as a function of 
UTC time in the lower left inset, (b) 1415 UTC with half hour of DC-8 flight track. 



Figure 1 shows a combination of two Geostationary Operational Environmental Satellite, GOES- 10 
and GOES-12, 10.8-pm IR brightness temperature images over the TC 4 domain at 1558 UTC, 5 August 
2007. Coordinated flights of the ER-2 over and the DC-8 in the clouds followed an east-west orientation 
in the Panama Bight at about 6.5°N, where an anvil cirrus was observed flowing from a deep 
convective core at about 79°W. One coordinated segment for the ER-2 and DC-8 over the anvil cirrus 
and deep convective core occurred from about 13:50 to 14:40 UTC (ER-2 flight) which is shown in 
Figure 2(a) together with the MASTER false color RGB image. The MASTER pseudo-RGB image is 
based on the band 5 (0.66 pm) reflectance, band 21 (2.16 pm) reflectance, and band 48 (11 pm) 
reversed-scale brightness temperatures. The ocean background is dark in the false color image. 
Optically thick ice clouds including the anvil and the upper part of the deep convective core are purple. 
Thin cirrus clouds are light blue. Note that only the DC-8 tracks (in red) coordinated with the ER-2 
track (in blue) are shown here. During the west-to-east leg, the DC-8 flight had a lag of about 7-8 min 
with respect to the ER-2, while along the east-to-west leg, the ER-2 lagged the DC-8 by about 6-7 min. 
Figure 2(b) shows the cloud vertical cross-section given by the CRS reflectivity along the ER-2 track. 
Cloud top and base detected by the CPL for the upper two ice cloud layers are also shown in Figure 2(b). 
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The vertical structure of the CRS reflectivity distinctly shows the outflow anvil cirrus from the 
deep convection. 

Figure 2. (a) MODIS/ASTER airborne simulator (MASTER) RGB images for flight 
tracks, 1350-1440 UTC, 5 August 2007 (blue: ER-2, red: DC-8, green: ER-2 during turn). 

(b) CRS reflectivity for the deep convective cloud system along the ER-2 flight track. 
CPL-detected cloud tops for the upper two layers are shown by the red and brown lines, 
respectively. CPL-detected cloud bases for the upper two layers are shown in black and 
gray lines, respectively. Gap corresponds to the ER-2 turn. 




14:00:00 14:12:00 14:24:00 14:36:00 

Time (UTC), 05 August 2007 


The DC-8 launched four dropsondes along the track shown in Figure 2: XI at 79.49°W, 6.50°N at 
13:54 UTC, X2 at 80.86°W, 6.19°N at 14:10 UTC, X3 at 79.68°W, 6.49°N at 14:23 UTC, and X4 at 
81.37°W, 6.18°N at 14:42 UTC. Surface radiosonde measurements launched at 15:05 UTC at Las 
Tablas, Panama (80.25°W, 7.75°N) from the Nittany Atmospheric Trailer and Integrated Validation 
Experiment (NATIVE) included high vertical resolution profiles of atmospheric temperature, 
humidity, and ozone mixing ratio [43], Figure 3 shows the mean profiles of atmospheric water vapor 
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and temperature from the four DC-8 dropsondes and the profdes of temperature and water vapor from 
the NATIVE ground-based sonde. The mean temperature profile is essentially the same as that from 
the NATIVE sonde. The mean water vapor profile from the dropsondes is only available below about 
10.3 km. Thus, above 10.3 km, the scaled NATIVE water vapor profile is used to complete the mean 
water vapor profile that is used for the present simulations together with the NATIVE temperature and 
ozone profiles. The dropsondes sampled two points, XI and X3, within the convective column and 
two, X2 and X4, outside of it under the anvil. Because the NATIVE sounding was taken in air lacking 
a deep convective cloud, a second sounding, launched at 23:20 UTC, 4 August 2007 from Alajuela, 
Costa Rica (9.98°N, 84.2 1°W) when a convective system was overhead, was used to better represent 
the convective cloud conditions that would be encountered in the convective cloud. The Alajuela 
sounding is significantly moister above 14 km than its NATIVE counterpart. 


Figure 3. Mean atmospheric water vapor and temperature profiles (with -IK shifting) 
from four DC-8 dropsondes (see text for details) along the track shown in Figure 2. 
Profiles from surface radiosondes launched at (a) 1505 UTC from the NATIVE site at Las 
Tablas, Panama (80.25°W, 7.75°N) and (b) at 23:20 UTC, 4 August 2007 from Alajuela, 
Costa Rica (9.98°N, 84.2 1°W). 
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3. Infrared Radiances over the Deep Convective Cloud System: Observations and Simulations 

3. 1. The S-HIS Derived IR Brightness Temperatures at the MODIS Bands 

Although the MASTER has bands similar to the MODIS complement, it has different spectral 
response functions (SRFs) and lacks a few channels. In particular, it does not have the MODIS water 
vapor bands at 6.7 and 7.3 pm. Here, the measurements from the S-HIS aboard the ER-2 are used to 
construct the MODIS bands at 3.7, 6.7, 7.3, 8.5, 10.8, and 12.0 pm. Figure 4 shows an example of an 
S-HIS spectrum and the Aqua MODIS and TC 3 4 MASTER band SRFs around 3.7, 6.7, 7.3, 8.5, 10.8, 
and 12.0 |tm. The S-HIS measurements and MASTER SRFs at 3.7, 8.5, 10.8, and 12.0 pm were used 
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to construct the pseudo-MASTER brightness temperatures, which were then compared to the 
corresponding observed MASTER brightness temperatures. The matched brightness temperatures, 
shown in Figure 5, agree well with each other with a standard deviation of ~1 K in their differences. 
Theoretically, if S-HIS and MASTER have the exact same observation geometry and FOVs, the 
comparisons should be essentially the same. The biases are likely due to calibration differences, while 
the variability in the differences is probably the result of the two instruments’ viewing slightly 
different scenes because of differences in their pixel sizes. At a cloud-top height of 15 km, the 
MASTER FOV is only ~12 m at nadir compared to 500 m for the S-HIS. Similarly, the S-HIS 
measurements and MODIS SRFs at 3.7, 6.7, 7.3, 8.5, 10.8, and 12.0 pm were used to simulate the 
MODIS-observed brightness temperatures at these bands (Figure 6). These pseudo-MODIS brightness 
temperatures are considered as MODIS observations in the present study. 

Figure 4. Example of S-HIS spectrum and Aqua MODIS and TC 4 MASTER band SRFs 

around 3.7, 6.7, 7.3, 8.5, 10.8, and 12.0 pm. (a) 3. 5-8.0 pm, (b) 8.0-12.8 pm. 



Figure 5. Brightness temperatures at (a) 3.75, (b) 8.66, (c) 11.3, and (d) 12.2 pm using the 
MASTER SRFs applied to S-HIS measurements and MASTER measurements. 
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Figure 5. Cont. 
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Time (UTC) Time (UTC) 


Figure 6. Equivalent MODIS brightness temperatures at (a) 3.79, 6.75, 7.33, and (b) 8.52, 
10.8, and 12.0 pm along the flight track constructed using MODIS SRFs applied to the 
S-HIS measurements. 



3.2. Optical Thickness Profiles from Combined CRS and CPL Measurements 

The vertical profile of the cloud water content can have significant influence on outgoing IR 
radiances, even for optically thick clouds. A thick convective cloud with a low ice water content at the 
top will produce a different IR temperature than one having a dense top (e.g., [44]). A multilayered 
ice-over-liquid cloud will yield a different IR temperature than a single-layered ice cloud having the 
same cumulative optical depth. Thus, it is important to use a realistic vertical profile of the water 
content or optical depth in the simulations. 

The CRS measurements are primarily used to obtain vertical cloud x profiles at a resolution of 1 5 
m. Three cloud phases, ice, liquid, and mixed phase, are considered here. The vertical profiles are 
separated into three layers, ice (top layer), mixed-phase (middle layer), and liquid (bottom layer) 
clouds. The layer with temperatures less than -20 °C is defined as ice cloud, the layer with 
temperature exceeding 0 °C is a liquid cloud. The layer between -20 °C to 0 °C is the mixed-phase 
cloud [45], In the bottom layer, the liquid cloud water content (L WC) is computed as follows, 

Z e = 57.54 LWC 5A1 (1) 

o /: 'i 

where ZIFC is in g-nf and Z e is in mm -m [46]. In the ice cloud, the ice water content (IWC) is 


Remote Sens. 2012 , 4 


3031 


IWC = 0.137 Z e 0 ' 643 (2) 

where IWC and Z e are in the same units as those for Z e -LWC [47], 

For the mixed-phase layer, LWC is assumed to decrease to 0 linearly with increasing altitude from 
the freezing level to the bottom of the ice cloud layer. Conversely, IWC is assumed to increase linearly 
from 0 g m -3 following Equation (2). On the basis of these assumptions, an ice fraction IceFra is defined 
to separate the two parts of CRS reflectivity, IceFra = ~(T - T 0 )/(T 0 - T Tra ), where T is atmospheric 
temperature in degree, To is the temperature at the freezing level, and T Tm is the transition temperature 
from mixed-phase to ice cloud. In the present study, 0 °C is used for To, so IceFra = T/T Tm . A phase 
transition temperature of -20 °C (altitude is ~8.3 km) is used. Therefore, the two parts of the CRS 
reflectivity due to cloud ice Z e we and liquid water Z e hq , are obtained from Zj ce = IceFra * Z e and 
Z t ! iq = (1 - IceFra) * Z e , respectively. IWC and LWC are then obtained by applying Z e hq and Z e hq to 
Equations (1) and (2), respectively. Because of difficulties in separating the surface radar return from 
the hydrometeors within the first 500 m above the surface, CRS reflectivity below 500 m is assumed to 
be that at the level of 500 m [48,49], 

While the CRS retrieval algorithm assumes the ice fraction varies as indicated above, the 
Rosemount icing probes detected no liquid water at temperatures up to -3 °C for this flight segment. 
Thus, to provide a more realistic variation of cloud water structure, a second profile was constructed 
assuming no liquid water for temperatures less than 0 °C. Thus, in the 0 to -20 °C transition zone used 
above, IceFra is set equal to 1.0 for a second set of calculations. Both results are used to determine the 
sensitivity of the calculations to the assumptions about supercooled liquid water in the clouds. 

The CPL measurements are used to compensate for the optically thin cirrus part missed by the CRS 
measurements [35], Figure 7 shows the vertical structure of the cloud system as seen from the CRS 
and CPL measurements and the 15-m resolution t profiles constructed from the CRS and CPL data for 
the pure ice layer. The total optical thicknesses, those from the CPL and those missed by the CRS 
along the flight track over the cloud system, are also shown in Figure 7. Visible x from the CRS shown 
in Figure 7(c) is derived from the empirical relationships, x = 0.065 -IWP 0M for cloud ice and 
x = 3.0/20.0 LWP for cloud liquid (assuming water clouds with D e = 20 pm and an extinction 
efficiency, Q e = 2.0), respectively [50], The water paths for each layer are the products of the water 
content and the 15-m thickness for each layer. The minimum detectable reflectance from the CRS is 
used to isolate the CPL-measured cloud areas missed by the CRS measurements: overlapping is 
denoted where the CRS gray region is visible under the CPL color region in Figure 7(a). The cloudy 
area missed by the CRS is shown as the gray region in Figure 7(b). Figure 8(a) shows examples of the 
constructed vertical profiles for 15-m layer ice and liquid cloud x. One profile is from an anvil cloud 
measured at 13:56:34 UTC, and the other is taken from a convective core measured at 14:06:10 UTC. 
The CPL cumulative ice cloud x missed by the CRS measurements along the flight path are detailed in 
Figure 8(b). The x value most frequently missed by the CRS but detected by the CPL is -0.20. Overall, 
the average x missed by the CRS is 0.28 ± 0.17, which is consistent with the values (0.15-0.45) 
deduced by [35], 



CPL Optical Thickness Altitude (km) Altitude (km) 
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Figure 7. (a) CRS-detected cloudiness and CPL-detected cirrus cloud layers, (b) vertical 
distributions of 15-m ice cloud CPL x and CRS x, (c) CPL-detected cirrus optical thickness 
and part missed by CRS measurements plotted with total optical thickness from CRS and 
CPL. The coincident DC8 tracks are shown in (a). Note the time difference of about 6-8 min 
between the DC-8 and ER2 tracks. 
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Figure 8. Cloud optical thickness summary, 5 August 2007. (a) Vertical, 15-m resolution x 
profiles at 13:56:34 UTC (anvil cloud, black) and 14:06:10 UTC (convective core, gray), 
(b) occurrences of CPL-derived cumulative ice cloud x in regions with no CRS detection 
(minimum detectable sensitivity). 
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3.3. Cloud Bulk Scattering Properties during the Flight Track 

Cloud bulk scattering properties are commonly used to simulate IR radiances emanating from 
clouds [51-55], Ice cloud bulk scattering properties are computed by averaging ice particle single 
scattering properties ( i.e ., single scattering albedo, absorption and scattering efficiency, asymmetric 
factor, and scattering phase function) over ice cloud habit and particle size distributions. 

Accurate measurements of ice cloud particle size and habit distributions are crucial for ice cloud 
bulk scattering properties. In this study, ice cloud particle size distributions from the 2D-S, CIP and 
PIP measurements are combined to construct the particle size distributions with a minimum size bin of 
5-15 pm. To minimize the effects of large-crystal shattering on number densities of small ice particles, 
the particle size distributions from the 2D-S for small ice particles are combined with those from the 
CIP and PIP for large ice particles. To be consistent with the CIP and PIP data that have been averaged 
over 5 s or about 1.0 km horizontal path, the 2D-S data are averaged in the same way. 

Figure 9 shows the combined ice particle size distributions over an anvil ice cloud and convective 
core from the 2D-S, CIP, and PIP measurements along the flight track over the deep convective cloud 
system (Figure 2). High concentrations of small particles are seen in both the anvil at 13:57:09 UTC 
(Figure 9(a)) and core at 14:02:07 UTC (Figure 9(b)). And both cloud types have a second peak 
concentrated at a maximum dimension D of 200-300 pm. All of the ice particle size distributions 
along the DC-8 flight track (Figure 7(a)) are derived and shown in Figure 10. The averaged particle 
size distribution is used for computing ice cloud bulk scattering properties. Besides the measured 
particle size distributions, the Gamma particle size distribution (e.g., [56]) is also used for investigating 
the effects of D e on the IR radiances from this cloud. The Gamma particle size distribution varies with 
its parameters, dispersion p and slope K. Those used by Hong [56] are adopted here as 

n(D) = NoD 2 exp(-4.87D/Do) 


( 3 ) 
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where No is the intercept and Do is the median of the distribution. 


Figure 9. Ice particle size distributions from combined 2D-S, CIP, and PIP measurements 
for (a) a case of cloud anvil (13:57:09 UTC), and (b) a case of convective core (14:02:07 
UTC) along the flight track. 



Figure 10. All ice particle size distributions from the 2D-S, CIP, and PIP measurements 
along the flight track. Mean particle size distribution shown in red. 



Also taken from the DC-8, the CPI data were used to resolve the ice cloud particle habits 
simultaneously with ice particle size distributions during TC 4 . Using an automatic crystal habit 
classification program, ice particle shapes measured by the CPI were placed into categories including 
sphere, column, plate, rosette, budding rosette, and irregular for ice particles with D < 50 pm and >50 pm. 
Figure 1 1 shows the ice particle habits detected by the CPI along the flight track. It was found that 
spherical and irregular particles comprise almost all of the habits when D < 50 pm, while the particles 
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with irregular and column habits are most common when D > 50 pm. The averaged habit mixing ratios 
along the flight track are used for computing ice cloud bulk scattering properties. When D < 50 pm, 
the habit distribution consists of 55% spheres, 37% irregular particles, and 8% columns. When D > 50 pm, 
the habit distribution consists of 73% irregular particles, 19% columns, 6% rosettes including rosettes 
and budding rosettes, 1% spheres, and 1% plates. Single-scattering properties of ice particles with 
various habits including column, hollow, bullet rosette, plate, aggregate, and droxtal have been 
investigated extensively (e.g., [53,55,57] and references therein). The single-scattering properties of 
droxtal, column, plate, rosette, and aggregate particles are used respectively for the sphere, column, 
plate, rosette, and irregular habits determined from the CPI data. 


Figure 11 . Ice particle habits detected from the CPI aboard the DC-8 along the ER-2 flight 
track during 5 August 2007 for ice particles with (a) D < 50 pm and (b) D > 50 pm. The 
averaged habit mixing ratios are also shown at the right. 




Time (UTC), 05 August 2007 


Following Yang et al., the single-scattering properties of ice particles are then averaged over the 
habit distributions and particle size distributions to obtain ice cloud bulk scattering properties, which 
are functions of wavelength, and D e , which is defined by: 
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where f[D) is the particle habit distribution and the summation of/(D) over N equals 1.0, i is the index 
for N ice particle habits at D, and V and A are the particle volume and projected area, respectively [53], 
The liquid cloud is assumed to be composed of spherical water droplets. An effective radius r e of 10 
pm with a gamma particle size distribution following Mishchenko and Travis [50] is used for the water 
cloud bulk scattering properties. For mixed-phase cloud portions (Figure 8 (a)), a water-phase mixing 
ratio y is defined by LWCKJWC + LWC) first, then the bulk scattering properties of mixed-phase 
clouds are derived by combining those of the ice and water clouds [58-61]. 

3.4. Effects of Optical Thickness and Particle Size on IR Radiances over Deep Convective Cloud 

Simulations and observations of the brightness temperatures at 3.7, 6.7, 7.3, 8.5, 10.8, and 12.0 pm 
(BT 03 , BToe, BToy, BTos, BTn, and BT 12 ) over the cloud system are compared to examine the effects of 
x and D e on the IR brightness temperatures. 

The correlated ^-distribution routines developed by Kratz [62] and Kratz and Rose [63] for the 
MODIS bands are used to compute the absorption optical depth for each layer in a clear-sky 
atmosphere. The atmospheric profiles of temperature and humidity (Figure 3), as well as ozone, are 
used in the correlated ^-distribution calculations for each band. The surface emissivity is assumed to be 
1.00 for the IR bands. The spectral brightness temperatures over the clouds are then computed with the 
discrete ordinates radiative transfer model [64] using the constructed cloud x profiles (Section 3.2), the 
computed bulk scattering properties, and various D e in the range of 10-150 pm. 

Figure 12 shows the simulated and observed BT values over the subject cloud system as a function 
of ice cloud x using the profile in Figure 3(b). Results are similar for most channels. Differences are 
discussed in Section 4. The mean particle size distribution ( D e = 95 pm) along the flight track is used 
for the simulations. Since the measured BT 03 values include reflected solar radiation while the 
simulations only include the emitted radiance, the measurements of BT 03 are omitted in Figure 12(a). 
The other simulated values generally agree with the corresponding observations. In particular, for 
optically thick ice clouds with x > 20 , the simulated values are essentially the same as the observations. 
For optically thin ice clouds (x < 1), the simulated BTos, BTn, and BT 12 are greater than the 
observations, while the simulated BT 06 , BT 07 , BTos, BTn, and BT 12 are less than the observations for 
ice clouds having 5 < x < 10 (anvil clouds). These differences are most likely due to the uncertainties 
in derived ice cloud x from the CRS (or/and CPL), which can be distinctly different because of the 
uncertainties in Z e -LWC relationship [44,65]. The underlying water clouds tend to decrease the BT 03 , 
BTos, BTn, and BT 12 values for non-opaque ice clouds, but have negligible effect on the water vapor 
channels’ brightness temperatures. The simulations and observations both reveal that all six BTs are 
sensitive to ice cloud optical depth for x up to 20. Furthermore, these sensitivities are essentially 
monotonic, i.e., these brightness temperatures monotonically decrease with increasing ice cloud optical 
thickness up to x = 20. When x is above 20, the brightness temperatures vary widely although a weak 
decreasing trend with increasing ice cloud x is apparent. 
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Figure 12. Simulated and observed above-cloud brightness temperatures at (a) 3.7, (b) 6.7, 
(c) 7.3, (d) 8.5, (e) 10.8, and (f) 12.0 pm as functions of ice cloud x using sounding in 
Figure 3(b), D e = 95 pm. 



0.1 1.0 10.0 100.0 0.1 1.0 10.0 100.0 
Optical Thickness Optical Thickness 


The sensitivity of the radiances to phase transition temperature is illustrated in Figure 13, where the 
3.7 and 6.7-pm simulations used the -20 °C transition temperature (red) and 0° transition (blue, same 
as blue in Figure 12). The results are similar for the other channels (not shown). Negligible change is 
seen in the radiances, but the entire set of radiances is shifted to the right. Thus, the simulated 
radiances correspond to greater optical depths if it is assumed that the hydrometeors are entirely in the 
ice phase for T < 0°, as indicated by the in situ data rather than assuming that the layer between 0° and 
-20 °C is a mixture of ice and water as formulated in Section 3.2 above. This result is encouraging 
because the phase of clouds in that temperature range is often uncertain. 

Figure 14 shows the simulated and observed BTDs among the 3.7, 6.7, 7.3, 8.5, 10.8, and 12.0-pm 
bands over the cloud system as a function of ice cloud x. As in Figure 12(a), the BTD(3.7-11) (BTD 
between 6.7 and 10.8-pm bands and others are similarly defined) measurements are not shown in 
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Figure 14(a). In general, the simulated values of BTD(6.7-11), BTD(8.5-11), BTD(11-12), and 
BTD(6.7-7.3) are consistent with their respective observations. They are differentially impacted by the 
presence of low clouds when the ice clouds are optically thin. The water vapor bands at 6.7 and 7.3 pm 
are affected by upper and middle tropospheric water vapor and above, respectively when the ice clouds 
are optically thin. This is clearly shown in Figure 14(f), which shows that the simulated BTD(6.7-7.3) 
values differ significantly from the observations when x < 3, but dramatically decrease as x increases. 
Both the simulated and observed BTDs show strong sensitivities to x of up to 10. Weak sensitivities to 
x up to 20 are also found for BTD(3.7-1 1), BTD(6.7-1 1), BTD(3.7-6.7) and BTD(6.7-7.3). 

Figure 13. Same as Figure 12, except only for (a) 3.7 and (b) 6.7 pm over the subject 

cloud as a function of ice cloud x using phase transition temperatures of 0 °C and -20 °C. 




Figure 14. Same as Figure 12, except for brightness temperature differences (BTD): 
(a) 3.7-11, (b) 6.7-11, (c) 8.5-11, (d) 11-12, (e) 3. 1-6.1, and(f) 6.1-1. 3 pm. 
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Figure 14. Cont. 
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Figure 15. Same as Figure 14, except dependence on effective particle size of ice cloud 
included: (a) 3.7-11, (b) 6.7-11, (c) 8.5-11, (d) 11-12, (e) 3.7-6.7, and (f) 6.1-13 pm. 
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Brightness temperatures for the six bands and their BTDs are commonly used to estimate ice cloud 
D e . As noted earlier, previous studies of infrared cloud retrievals have been limited mainly to non-opaque 
ice clouds having x < ~5. Here, the effects of ice particle size D e on the BTDs and the potential of 
estimating D e for optically thick ice clouds are investigated by performing simulations with D e in the 
range of 10-150 pm. Figure 15 shows simulated BTDs for various 3.7, 6.7, 7.3, 8.5, 10.8, and 12.0-pm 
pairs over the subject cloud system as functions of x and D e . For these simulations, only ice clouds are 
involved. The observed BTDs, except for those for the 3.7-pm band with its reflected solar radiances, 
are also shown together with these simulations. Consistent with previous studies, BTD(3.7-11), 
BTD(8.5-11), BTD(11-12), and BTD(3.7-6.7) show strong sensitivity to ice cloud D e when the ice 
clouds are non-opaque [10,14,19,66], With x increasing above ~3, the differences in BTD(3.7-11), 
BTD(8.5-11), BTD(11-12), and BTD(3.7-6.7) as a function of D e decrease eventually. BTD(6.7-11) is 
essentially insensitive to D e , which indicates that it could provide a good estimation of r, since it 
increases asymptotically toward zero at r> 20 in the present simulations. BTD(6.7-7.3) only shows 
differences for D e < 30 pm. When the opaque ice cloud x increases, the sensitivity of these BTDs to D e 
decreases dramatically. BTD(8.5-1 1) and BTD(1 1-12) are only sensitive to small D e while BTD(3.7-1 1) 
and BTD(3.7-6.7) still show weak sensitivity to large values of D e . 

Figure 16. Brightness temperature differences for select pairs of 3.7, 6.7, 8.5, 10.8, and 
12.0-pm radiances over the tropical deep convective clouds (20 < x < 100) as a function of 
effective particle size D e . RSD is relative standard deviation. 



Although BTD(3.7-11), BTD(8.5-11), BTD(11-12), and BTD(3.7-6.7) show significant variability 
over the deep convective clouds (x > 20), the general feature of decreasing BTDs with increasing D e is 
still apparent. In order to investigate the potential for estimating D e , the BTDs for each D e are averaged 
for ice clouds having x in the range of 20-100. Figure 16 shows the BTD(3.7-11), BTD(8.5-11), 
BTD(11-12), and BTD(3.7-6.7) over the tropical deep convective clouds (20 < x < 100) as a function 
of D e . It is evident that BTD(3.7-11) and BTD(3.7-6.7) are sensitive to D e , especially for D e > 50. 
BTD(8.5-1 1) is only sensitive to D e values smaller than 50-70 pm and BTD(1 1-12) is only sensitive to 
the smallest D e values. This feature reveals the potential for estimating small D e using those two BTDs 
for tropical deep convective clouds during nighttime. The dependencies of BTD(3.7-11) and 
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BTD(3.7-6.7) are likely too small to be useful for a retrieval of D e for x > 20 for two reasons. At the 
low temperatures associated with these thick clouds, the noise at ~3.7 pm can be as high as 3.5 K [67] 
and the uncertainty in the water vapor in and around the cloud top can cause errors of ~1 K (shown 
later), which cover the range of the signal in Figure 16. 

4. Discussion 

Simulations of cloud radiance fields necessarily require many assumptions about the cloud, 
humidity, and temperature horizontal and vertical variations. Differences between the observations and 
simulations in Figure 14(b,f) are likely due to those assumptions. As noted above, the differences for 
x < 6 are probably due to the presence of low or midlevel clouds below thin cirrus clouds. For larger 
optical depths computed from the CPL and CRS data, the observed water vapor channel BTDs do not 
appear to reach a distinct asymptote around x = 8 as seen in the simulations. Rather, the asymptote is 
not reached until x = 20 or so. Thus, the vertical variations in D e , IWC, temperature, and humidity are 
likely different than assumed and/or higher layer resolution is needed in the calculations. Given the 
complexities in the actual structure, however, the results are remarkably similar overall. 

4.1. Sensitivity of BTDs to Uncertainties in Atmospheric Profiles 

Nevertheless, it is important to understand how errors in the atmospheric profiles affect the modeled 
BTDs. The results in Figures 12-15 were all based on the sounding in Figure 3(b), which was 
considerably moister in the upper troposphere than the profde in Figure 3(a). Such differences may be 
representative of errors in numerical weather analyses where tropical convection may be properly 
diagnosed, but offset spatially from where it is actually occurring. To determine how the temperature 
and moisture profiles affect the radiances, the calculations were performed using the composite smooth 
sounding in Figure 3(a). The differences in BT and BTD for the two sets of calculations are very small, 
except for the water vapor channels and their BTDs. Figure 17 shows results for two water vapor 
channel pairs computed using the two profiles in Figure 3. In all cases, BTD(6.7-1 1) and BTD(6.7-7.3) 
are greater using the moister Alajuela sounding. The discrepancies in the BTDs due to the profile 
difference are mostly between 0.5 and 1 K for x > 10. 

To gain further insight into the sensitivity of the radiances to such errors, four different experiments 
were conducted. The BTDs of interest were recomputed by systematically varying the upper 
tropospheric humidity (UTH), the upper tropospheric temperature profile, the cloud top height, and the 
ice water content profile for the observed cloud system. Results for BTD(3.7-11), BTD(8.5-11), and 
BTD(11-12) show negligible changes for the scaling in each case. The results for BTD(3.7-6.7) are 
comparable to those for BTD(6.7-1 1). Thus, only BTD(6.7-1 1) and BTD(6.7-7.3) are shown. 

Figure 18 shows changes in BTD(6.7-11) (left) and BTD(6.7-7.3) (right) due to systematic 
variations in water vapor (WV) above 7 (top) and 10 km (bottom) along the aircraft flight path. 
Because the ice cloud optical thickness varies along the track (Figure 7(c)), the changes vary with 
x (replotted in Figure 18(d)). Each curve indicates the differences between the BTDs computed using 
the nominal relative humidity scaled by 0.5, 0.75, 1.25, and 1.5 and those computed using the nominal 
humidity profile in Figure 3(b). Decreasing (increasing) the UTH increases (lowers) the BTDs for both 
spectral pairs. The impact is greatest for thin clouds, particularly for BTD(6.7-1 1). For clouds having 
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x > 10, the change in BTD < 1 K, unless a thin cloud lay above the thicker cloud. Overestimation of 
the UTH for the thin clouds can help explain why the BTDs in Figure 14(b,f) are low relative to the 
observations. Obviously, the vertical profile of humidity will be considerably different in a column 
completely filled with clouds than one having clouds in only one layer of the column, so UTH errors 
of 50% or greater in the vicinity of a cumulonimbus would not be surprising. The impact of the 
humidity on the BTDs for other band pairs was negligible, except for BTD(3.7-6.7), which yielded 
results similar, though reversed in sign, to those for BTD(6.7-1 1). 

Figure 17. Same as Figure 14, except only for (a) BTD(6.7-11) and (b) BTD(6.7-7.3) 

computed using soundings in Figure 3(a) (Profile 1) and in Figure 3(b) (Profile 2). 




Figure 18. Changes in BTD along flight track relative to BTDs computed using sounding 
in Figure 3(b) for various scalings of the sounding’s water vapor (WV). (a) BTD(6.7-11) 
and (b) BTD(6.7-7.3) for scaled WV above 7 km, (c) BTD(6.7-11) and (d) BTD(6.7-7.3) 
for scaled WV above 10 km. 
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Similar scaling was applied to the IWC above 10 and 12 km and used to compute the BTDs. 
Changing the IWC effectively alters the optical thickness of the cloud layer, so it should have a 
significant impact on thin clouds. Figure 19 presents the results for BTD(6.7-11) and BTD(6.7-7.3). 
The impact is less than 1.1 K for x > 20 and smaller than 2 K for 10 < x < 20. For those cases, an error 
in IWC of 100% is needed to raise the difference above 1 K. For x < 10, halving the IWC drops the 
BTD(6.7-11) by as much as 15 K, while doubling it increases the BTD by as much as (Figure 19(a)). 
The impacts are much smaller when only part of the cloud is altered (Figure 19(c)) and for BTD(6.7-7.3) 
as seen in Figure 19(b,d). Biases in the IWC profile could also contribute to the relatively large 
differences between the computed and observed BTDs seen in Figure 14(b,f). For optically thick 
clouds, the variation in IWC had an impact less than 0.1 K on BTD(1 1-12) (not shown). BTD(3.7-6.7) 
changes due to the IWC variations were slightly larger than those for BTD(6.7-1 1). 

Figure 19. Same as Figure 18, except for various scaling of the ice water content for 
(a) BTD(6.7-1 1) and (b) BTD(6.7-7.3) above 10 km and for (c) BTD(6.7-11) and 
(d) BTD(6.7-7.3) above 12 km. 
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Altering the temperature profiles between 10 and 15 km and between 12 and 15 km by 2 or 4 K 
produced changes in BTD(6.7-11) and BTD(6.7-7.3) of less than 0.4 K for 10 < x < 20 and less than 
0.2 K for x > 20 (not shown). Moving the cloud higher by as much as 3 km produced an increase in 
both BTD(6.7-11) and BTD(6.7-7.3) with values ranging from 0 K for optically thin clouds to 
increases of up to 4 K for 5 < x < 20 and down again to ~ 1 K for x > 20 (not shown). The increase in 
BTD(6.7-11) is consistent with a stratospheric penetration by convective clouds [68]. Lowering the 
cloud by 3 km decreases the thick cloud (x > 20) 6.7-pm BTDs by as much as 6 K and by up to 15 K 
for mid-range optical depths. This decrease results from increased WV above the cloud. 
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The sensitivities illustrated above contribute to uncertainties in the radiance calculations relative to 
the observations, particularly for the unsampled portions of the cloud contributing to the observed 
radiances. Differences between the various sensor fields of view also introduce some variability in the 
observed radiances relative to the calculations. For the most part, it appears that variations in D e near 
cloud top could explain much of the thick cloud BTD variability in Figure 15(c,d), but not for the 
water vapor channel BTDs in Figure 15(b,f). For those channel pairs, a combination of sensor 
calibration errors, sensor viewing angle differences, and uncertainties in the input data are likely to be 
responsible for 1-2 K variability in the observations that is not seen in the calculations. 

The above results suggest that there is useful information about thick ice cloud optical depth in the 
BTDs between several channels commonly found on operational meteorological satellites. In 
particular, the 6.7-pm water vapor channel, when differenced with the 11 -pm channel, appears to be 
most sensitive to x. The sensitivities to errors in the soundings, cloud temperature, and IWC studied 
here are typically on the order of 1 K for optically thick clouds. Thus, while the sensitivities of 
BTD(6.7-11) to changes in thick ice cloud x indicate that methods to retrieve x at night can be 
developed, their results will likely be subject to uncertainties greater than those encountered using 
solar reflectance methods. 

4.2. Satellite Examples of T Dependency on BTD(6. 7-11) 

While it is beyond the scope of this paper to completely evaluate the typicality of the above results, 
it is useful to determine if the dependencies seen here are, at least, qualitatively evident in satellite 
imagery. Figure 20 shows imagery over Central America from GOES- 12 taken at the same time as the 
example presented above. The pseudocolor RGB (red: BTn, scale reversed; green: BTD(3.9-11); and 
blue: 0.65-pm reflectance) and BTD(6.7-11) images are shown in Figures 20a and 20d, respectively. 
Figure 20(b,c) depicts the effective cloud height Z c and x, respectively, derived using the methods of 
Minnis et al. [10]. Generally, in these images, wherever Z c > 6 km and x > 4, BTD(6.7-1 1) > -20 K or 
so. While there is not perfect correlation, some of the large optical depths are caused by thin cirrus 
over thick stratus, it is clear that the gradients in BTD(6.7-11) correspond, to some extent, to the 
gradients in x for the thicker portions of the high clouds. 

At night, similar gradients in BTD(6.7-1 1) are evident around the convective cells in the 1045 UTC 
GOES- 12 imagery for 5 August 2007 (Figure 21). The correspondence between the highest clouds and 
the BTD(6.7-11) gradients around the cores is as strong in Figure 21(b,d) as those seen in Figure 20, 
suggesting similar variations in x. The imagery also highlights some of the difficulties that will be 
encountered when trying to use such spectral signals to retrieve x. The midlevel (5-7 km in Figure 21(b)) 
clouds over the Andes and even some of the marine stratus clouds in the lower left part of the image 
are also accompanied by some variations in BTD(6.7-11) that approach -20 K. With even drier 
atmospheres over mid- and lower-level clouds in other areas, BTD(6.7-1 1) will be even closer to 0 K. 
Thus, other information such as cloud height (e.g., Figure 21(b)) and other channel BTDs, as well as 
accurate temperature and water vapor profiles would be needed to unravel the BTD(6.7-1 1) signal. 
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Figure 20. GOES-12 imagery and retrieved cloud properties, 1415 UTC, 5 August 2007. 
(a) RGB, (b) effective cloud height, (c) optical depth, (d) BTD(6.7-1 1). 



There are considerable variations in BTD(3.9-11) in Figure 21(c) with the largest values (>30 K) 
occurring in the anvil-generated cirrus. The 3.7 and 3.9-pm channels are similar enough to make some 
inferences from Figure 15(a). The greatest values likely correspond to smaller values of D e for 2 < x < 8. 

As the optical depth increases toward the core, BTD(3.9-11) drops precipitously to <10 K and 
eventually becomes noisy over the core because the accuracy of BT 0 3 is greatly diminished at the 
lowest temperatures (<205 K). The gradient in BTD(3.7-11) over a convective core can be seen more 
clearly for the developing storms just southwest of Hispaniola. There, BTD(6.7-11) is close to -1 K, 
which would correspond to 30 < x < 100 (Figure 15(b)) and the smallest value of BTD(3.7-1 1) is ~5 K, 
which would correspond to small particles at the cloud top (Figure 21(b)). While the earlier discussion 
regarding the small gradient in BTD(3.7-1 1) with increasing D e would be overwhelmed by noise in the 
detector at low temperatures, the satellite observations suggest that D e could be retrieved at night for 
many of the thicker ice clouds, but not for the coldest clouds. More likely, such a retrieval would be 
limited to clouds with x < 20. Further research is needed to clarify the range of optical depths and 
temperatures for which such retrievals would be useful. 

The case presented in the previous section used a tropical convective system and may have 
characteristics unlike high clouds in mid-latitudes and polar regions. Figure 22 shows imagery from 
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the Meteosat-9 Spinning Enhanced Visible Infrared Imager (SEVIRI) taken at 1200 UTC, 1 June 2012 
over Europe. The RGB image (Figure 22(a)) shows baroclinic systems over the eastern and western 
sides of the domain with stratus clouds between and cirrus clouds over Italy and North Africa. Cloud 
optical depths retrieved with the methods of Minnis et al. [10] exceed 100 over parts of the eastern 
system and are greater than 20-40 over portions of the western high clouds. Although the 
correspondence of BTD(6.7-11) in Figure 22(c) with x is less obvious than seen in Figure 20, some is 
still present. This diminished correspondence could be due to the considerable amount of multi-layered 
clouds occurring in the baroclinic systems and the higher viewing zenith angles for the domain. The 
variations in x derived using the 0.65-pm reflectance are due to the total column cloud optical depth, 
while x for the upper layer clouds should primarily be responsible for variations in BTD(6.7-11). 
While this difference in sensitivity could raise some obstacles to retrieving x at night, it could also 
provide a means for detecting multilayered clouds and, perhaps, retrieving their properties using an 
approach similar to that of Chang et al. [69]. 


Figure 21. GOES- 12 imagery and retrieved cloud properties, 1045 UTC, 5 August 2007. 
(a) RGB, (b) effective cloud height, (c) BTD(3.9-1 1), (d) BTD(6.7-1 1). 
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Figure 22. Meteosat-9 SEVIRI imagery and retrieved cloud properties over Europe, 1200 
UTC, 1 June 2012. (a) pseudocolor RGB image, (b) retrieved cloud optical depth, 
(c) BTD(6.7-1 1), (d) BTD(8.5-1 1), (e) BTD(11-12), and (f) multilayer cloud index: thin 
cirrus over low water: magenta, likely thicker cirrus over low water: yellow, likely 
contiguous thick ice over thick water: brown. 



The information in the BTD(8.5-1 1) and BTD(1 1-12) images in Figure 22(d) and 22(e), respectively, 
could also aid in the determination of multilayered cloud conditions and the utility of the BTD(6.7-1 1) 
for thick cloud T retrievals. For example, the large values of BTD(8.5-1 1) and BTD(11-12) over Italy 
and North Africa indicate that the high clouds have T < ~4 and small values of D e . The near-zero 
values for both channel differences over much of the frontal cloudiness suggest that the ice clouds are 
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opaque and the retrieval of x using BTD(6.7-11) would apply only to ice clouds. If it is significantly 
different from the visible channel x retrieval, then it could be concluded that a multilayered or 
contiguous, very thick cloud was present. Figure 22(f) shows the results of applying the 1 1 and 13.4-pm 
channel method of Chang et al. to the same image [69], Most of the clouds classified as thick multilayer 
or thick high clouds by the Chang et al. technique correspond to areas where BTD(6.7-1 1) > -15 K [69]. 
Thus, it appears that even in the absence of the 13.4-pm channel, it may be possible to obtain valuable 
multilayered cloud information using the 6.7 and 11 -pm channels common to all current and most 
older geostationary satellites. 

5. Conclusions 

With the aid of TC 4 aircraft measurements of radiances over a tropical deep convective cloud 
system and simulations of radiances using in situ measurements of cloud properties, this paper has 
investigated the sensitivities of brightness temperatures for these bands and their differences to ice 
cloud x and D e . Measurements from two aircraft and surface-based radiosondes were used in a detailed 
radiative transfer model to simulate radiances at various MODIS channels that would emanate from a 
deep convective cloud complex. The simulated radiances were compared with radiance measurements 
taken at nadir from two instruments on the ER-2 and adjusted to the MODIS filter functions. 

The simulated BTs are generally consistent with the corresponding observations, particularly, for 
optically thick ice clouds with x > 20 (i.e., deep convective clouds). The IR radiances emitted from 
deep convective clouds are essentially the result of contributions from the cloud and from atmospheric 
gaseous absorption and emission above the clouds (upper troposphere), which can be neglected for the 
window bands, 3.7, 8.5, 10.8 and 12.0 pm. Both the simulations and observations of these IR 
brightness temperatures monotonically decrease with increasing ice cloud x up to ~20. Although, in 
some instances, significant differences are found between the absolute values of the simulations and 
observations, a weak decreasing trend is found with increasing ice cloud x when x is above 20. Similar 
to the brightness temperature comparisons, simulated BTDs among these bands generally agree with 
the observations within the range of uncertainties in the input data. Both simulated and observed 
BTD(3.7-11), BTD(6.7-11), BTD(3.7-6.7), and BTD(6.7-7.3) are sensitive to ice cloud x up to ~20. 
This dependency is accompanied by uncertainties due to sensitivities of the BTDs to errors in the 
atmospheric profiles of temperature, humidity, ice water content, and cloud top height that can be as 
large or larger than the variation of BTD with x. 

The sensitivity of BTDs to D e was also investigated. BTD(6.7-11) has no sensitivity to D e . This 
feature provides an advantage for estimating ice cloud x. Consistent with previous studies on 
estimating ice cloud D e using IR BTDs, BTD(3.7-11), BTD(8.5-11), BTD(11-12), and BTD(3.7-6.7) 
show strong sensitivity to D e of non-opaque ice clouds. With increasing ice cloud x, the sensitivity 
weakens. Over deep convective clouds (x > 20), both BTD(3.7-1 1) and BTD(3.7-6.7) show sensitivity 
to D e while BTD(8.5-1 1) is only sensitive to small D e that are less than 50-70 pm and BTD(1 1-12) is 
sensitive to even smaller D e . Unfortunately, the D e signal is likely to be overwhelmed by the noise in 
some channels at cold temperatures and uncertainties in the vertical profiles of temperature, humidity, 
and ice water content. Additional research is needed to define the range of temperatures and optical 
depths for which D e could be retrieved using the spectral bands considered here. 
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The simulated results are also generally consistent with satellite observations taken at similar 
wavelengths and indicate that some combination of these channels, particularly the 3.7, 6.7, and 1 1-pm 
channels, could be used for retrieving the optical depths of ice clouds that exceed the typical opaque 
limit of 6 for cloud optical depth. While many challenges and obstacles, such as the size of the 
uncertainties in input data, would need to be overcome, it may be possible to use such information to 
provide better quantitative monitoring of ice clouds over the entire diurnal cycle. Such a capability 
would aid both climate and weather studies because clouds are an essential component of the radiation 
and water budgets of the atmosphere both day and night. 
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